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Abstract 



We present a self-contained development of Gauss' Arithmetic-Geometric Mean 
' (AGM) and the work of A. E. Ingham who obtained rigorous error bounds for the 

■ AGM's approximations to the period of a simple pendulum. 

1 Introduction 

One of the most celebrated problems in the classical dynamics of particles is the computation 
of the period of the simple pendulum. The nonlinear differential equation which models the 
; pendulum's motion appears in numerous physical problems and the exact formula for the 

period of a single oscillation is given by a complete elliptic integral of the first kind. But, 
^ ■ in 1834, Joseph Liouville [9] proved a justly famous theorem which implies that such an 

! integral cannot be evaluated by any finite combination of elementary functions. Therefore, 

the calculation of the period must be carried out by suitable approximative formulas. This 
had been recognized long before, and in 1747, Daniel Bernoulli [12] published the first 
such approximation. Since then an enormous literature has arisen around the problem of 
finding a good approximation to the period and research continues unabated to this very 
r£ i day. 

The authors of these approximations show great dexterity and ingenuity in their deriva- 
tions and use a variety of techniques to obtain them. However, virtually NONE of them 
offers a rigorous error analysis. That is to say, inequalities on the upper bound for the 
error, which shows how good the approximation is, and on the lower bound for the error, 
which shows how bad the approximation is. (See Thurston [H]). Most of the authors do 
include numerical studies of the accuracy of their approximations and some even include a 
few order-of-magnitude asymptotics. But those with rigorous error bounds are few and far 
between. 

Recent interest has concentrated in Gauss' Arithmetic-Geometric Mean (AGM) algo- 
rithm [5] because of its high rate of convergence. In 2008, Claudio G. Carvalhaes and 
Patrick Suppes [2] published a very interesting and detailed presentation of the AGM and 
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its application to the approximation of the period. They also presented an elegant interpre- 
tation of the AGM recurrence formula as a method of renormalizing the pendulum in the 
sense that it replaces the original pendulum with another one with the same period, but 
longer length and smaller amplitude. This interpretation was already known to Greenhill 
[6] in the late 1800 's, but has been woefully neglected till recently. However, their paper, 
too, fails to offer any rigorous error analysis, although the numerical studies of the error are 
extremely interesting and merit study. 

It is unfortunate that none of the authors cites the marvelous investigations of the great 
British number theorist A. E. Ingham which L. A. Pars describes in his monumental 
665-page standard work which was published almost 50 years ago in 1965. Ingham 
not only obtains the formulas of Carvalhaes and Suppes but also obtains rigorous error 
estimates, both in excess and in defect. It is beyond question that Ingham's work deserves 
to be better known. 

So, our paper is organized as follows. To make it as self-contained as possible, we develop 
ab initio the theory of the AGM including Gauss' original proof that it converges to the 
complete elliptic integral of the first kind. It is difficult to find this anywhere, today, since 
the clever method of D. J. Newman [TU] has now become fashionable. Then we apply 
the AGM to the case of the simple pendulum and we slightly alter the results and proofs 
of Ingham so as to obtain a complete error analysis. Finally we compare our analytical 
error bounds with the numerical studies of Carvalhaes and Suppes and show that they 
virtually coincide (as they should!). But, our analysis explains the why of their numerical 
results. 

2 The AGM 

We use the Stockholm lectures of Vladimir Tkachev [U] in our treatment of the AGM. 

GAUSS' only published account of his arithmetic- geometric mean (AGM) algorithm ap- 
pears in a paper |5j on secular variations published in 1818. Yet, his unpublished papers [5] 
show that he was well aware of its major properties by 1799. 

Definition 2.1. Let a ^ and b ^ be two numbers such that a ^ b ^ and define the 
numbers ao and bo by 

a := a, b := b; (2.1) 
then for n — 0, 1, 2, . . . define the sequences {a n } and {b n } by 

fln+l := ~^~2 — ~> := y^n^n- (2-2) 

Note that each a n+ i is the arithmetic mean of the previous a n and b n , while each b n+ \ is 
the geometric mean of those same two numbers. 

Definition 2.2. One says that the sequences {a n } and {b n } in (12. ip and (12. 2 p define the 
arithmetic-geometric mean algorithm, which we abbreviate as AGM. 
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Now we collect some of the basic properties of the AGM. 

Proposition 2.3. The following properties of the AGM are valid. 

1. The a n 's decrease, the b n 's increase and every a n is bigger than every b m . More pre- 
cisely, 

ao ^ a l ^ • • • ^ a n ^ a n+1 ^ • • • ^ b n+1 ^ b n ^ • • • ^ b 1 ^ b . (2.3) 

2. 

0^a n -b n ^ . (2.4) 

3. The limits 

A := lim a n and B := lim b n (2-5) 
both exist and they are equal, 

A = B. (2.6) 

Proof. Of Q: 

Since the square of a real number is always non-negative, it follows that for n = 0,1,2, ... , 
(-v/a^ — v^n) 2 ^ and that there is strict inequality unless a n = b n , whence we conclude 
that the following inequality is valid, 

a n + K 



> Va~J^. (2.7) 

Of course, (12.71) is the famous arithmetic- geometric mean inequality for two numbers. Ap- 
plying it to a n+ i and b n+ i we obtain 

a n +i > b n+ i. (2.8) 
Thus, from a n+ i ^ b Tl+ \ and a n ^ b n we obtain 

a n > ±±± =: a n+1 > b n+1 := yfeK > b n , (2.9) 

which is fTOj) . 

of m-- 

From 6 n+ i ^ b n we conclude 

a n+l — o n+1 ^ a n+1 — o n = o n = 

and (I2.4p follows by induction. 
Of ([23D and (1231) : 

By ( 12 .3p the sequence {a n } decreases monotonically and is bounded from below by b , and 
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so A exists. By (12. 3 j) the sequence {6 n } increases monotonically and is bounded from above 
by a , and so B exists. 

Finally, letting n tend to infinity in ( 12. 4 p and using (12. 5p . we obtain 

A- B < 

and by the "squeeze" theorem, we conclude A = B. □ 
Now the following definition makes sense. 

Definition 2.4. We define the arithmetic-geometric mean, M(a, b) = /x of the numbers 
a and to be the common limit 

M(a, b) = /I := A : = lim a n = B := lim b n (2.10) 

n— >oo n— >oo 

of the AGM as applied to the numbers a and b. 

Proposition 2.5. The geometric mean b n is a closer approximation to \i than a n ; more 
precisely 

< - - bn < 1. (2.11) 
a« - A* 

Proof. We observe 

a n + b n . 

/i < a n+1 = — - — 2/i < a n + 6 n /j, - b n <a n - /i. 

Since < /i — b n < a n — fi, we can divide by a n — /i to complete the proof. □ 
3 Gauss' theorem on elliptic integrals 

The following theorem gives a hint of the depth of the mathematics involved in the AGM. It 
is the only theorem Gauss published on the algorithm and appears in the 1818 paper already 
cited [5]. But, it seems that he already had a proof in 1799 [3]: 

Theorem 3.1. Let a and b be positive real numbers. Then 
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f* d<t> 


M(a, b)~ nj 


o \J a 2 cos 2 <p + b 2 sin 2 (ft 



(3.1) 



The integral (13. ip is a complete elliptic integral of the first kind and, as we have already 
seen [9], cannot be evaluated in finite terms with elementary functions. In the next section 
we will see its relationship to the simple pendulum. 

Before we enter into the details of Gauss' proof, we introduce some notation and separate 
out the fundamental technical step. 
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Let 



I{a, b) :-- 



o a/ a? cos 2 + b 2 si 



sin 



Then, we have to prove that 

/(a, 6) = /(ai, &i) = I(a 2 , b 2 ) = I(a 3 , 6 3 ) 
since we can then conclude that 



7T 



I(a, b) = lim I(a n ,b n ) = I(/x,/x) = — 

n— >oo ZlJL 



(3.2) 



(3.3) 



(3.4) 



which, after multiplying by — , is precisely (13. ip . 

7T 

In order to conclude 



lim I(a n , b n ) = /( lim a n , lim 6„) = /(/x, /x) 



(3.5) 



we have to prove that we can interchange the limit and the integral signs. For this, it is 
sufficient to prove: 



Proposition 3.2. The sequence 



formly to —. 

H 



a 2 cos 2 4> + f> 2 sin 



, n = 0, 1, 2, . . . , converges uni- 



Proof. That means given any e > we must prove there exists positive number iV(e), which 
is independent of the variable (/), such that the following implication is true: 



n > N(e) 



1 



< e. 



a/ a 2 cos 2 <p + b\ sin 2 <fi ^ 
However, the identity cos 2 <fi + sin 2 <fi = 1 as well as the inequalities (I2.3P and 



(3.6) 



b n ^ v a 2 cos 2 + 6 2 sin ^ a T , 



show us that 



a n - K) = b n - a n < b n - /x < J a 2 cos 2 + 6 2 sin - /i < a n - /x < a n - b n , 



that is, 



' a 2 cos 2 + fr 2 sin 2 — /x 
where we applied (12 ,4p in the last inequality. Now 



< a n - b n < 



a — b 
2 n 



(3.7) 



by (13. 7p and (12.31) . For the implication (13. 6p to be true, it is sufficient that the following 
inequality be true: 

a-b l a-b ln(g) 

2™ b 2 <€ b 2 e ^ n> i n2 ' 

that is, the choice 

N(e) := (3.8) 



proves the truth of the implication ( I3.6p . and that, therefore, we can interchange the limit 
and integral signs in (13. 5ft . □ 

Gauss' original proof is based on the following change of variable in the integral I(a,b): 
we introduce a new variable, 0' instead of by the formula: 



2a sin 0' 

sin i 



a + b + (a — b) sin 2 0' 



(3.9) 



Proposition 3.3. Under the mapping (13.91) the interval ^ 0' ^ | corresponds bijectively 
to the interval ^ ^ ~ . 

Proof. Define the function 

/(*):= 2 ^ — - . (3.10) 

J w a + 6 + (a - o)t 2 v y 

Then 

{a + o + (a — 6)^}^ {a + b + (a — o)t } 
which proves that f(t) is increasing on [0, 1]. Moreover, 

/(0) = 0, /(1) = 1, 
which shows that /(£) maps [0, 1] bijectively onto itself. This completes the proof. □ 

Proof of Gauss' theorem on elliptic integrals. Gauss, himself j5], first states Proposition 13.31 
Then he blithely asserts u Evolutione autem rite facta, invenitur esse" which translates to 
" After the development has been made correctly, it will be seen (that). . . 

(3.12) 



a/ a 2 cos 2 + b 2 sin 2 \J a\ cos 2 0' + b\ sin 2 0' 

(We have changed Gauss' notation: he writes m, n, m' , n', T, X" in place of our a, b, ai,b\, 0, 0', 
respectively.) This, of course, is the step 

I(a, b) = I (a x , h) (3.13) 
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in the notation (13. 2p . 

We will show how the development is "made correctly" (!). 

Claim 1: 

2 cos 0' a/ a\ cos 2 0' + b\ sin 2 M 
a + b + (a — b) sin 2 0' 



cos0 = j\ ~ 2 // • (3-14) 



Proof. By and 1122]) . 



cos 2 0=1 — sin 2 - 



4a 2 sin 2 0' 



{(a + 6) + (a - 6) sin 2 0'} 2 
(a + 6) 2 + 2 (a 2 - b 2 ) sin 2 0' + (a - fe) 2 sin 4 0' - 4a 2 sin 2 0' 
~ {{a + b) + (a -6) sin 2 0'} 2 

4a 2 - 4(2a 2 - fr 2 ) sin 2 0' + 4(a 2 - of) sin 4 0' 
~~ {(a + b) + (a-b) sin 2 0'} 2 

_ 4a 2 cos 4 0' + 46? sin 2 0' cos 2 0' 
= {(a + 6) + (a -6) sin 2 0'} 2 

and factoring out 4 cos 2 0' and taking the square root of both sides gives us (13.1 4ft . □ 

Claim 2: 



(a + b) + (a- b) sin 2 A ' 
(a + b) — (a — 6) sin 2 0' 



a 2 cos 2 + 6 2 sin 2 = a T.,]iT~ Z* T, ■ (3-15) 



Proo/. By fl3~T41) . (ESj) and CT . we obtain 



9 9 , , 9 9 , 9 f 2 cos 'a/o? cos 2 0' + 6? sin 2 0' 1 4a 2 6 2 sin 2 0' 

a 2 cos 2 + b 2 sin 2 = a 2 <^ r v 1 y . | ^ + ; . 2 — 

i a + b+ (a - 6) sm 2 0' J {(a + 6) + (a - b) sin 2 0'} 2 

4a 2 cos 2 0'(a? cos 2 0' + b\ sin 2 0') + 4a 2 b 2 sin 2 0' 

{(a + 6) + (a-b) sin 2 0'}' 



I \ U ~T U) -r- | (i — i; / mil (V J" 2 

_ 2 a?(l - sin 2 0') 2 + b\ sin 2 0'(1 - sin 2 0') + 6 2 sin 2 0' 



4 a " 

{(a + 6) + (a-b) sin 2 0'} 2 

2 (a + 6) 2 (1 - sin 2 0') 2 + 4a6sin 2 0'(1 - sin 2 0') + (a - b) 2 sin 4 0' 

a {(a + 6) + (a - 6) sin 2 0'} 2 

2 (a + 6) 2 - 2(a - 6) (a + 6) sin 2 0' + 46 2 sin 2 0' 



{(a + 6) + (a-6) sin 2 0'} 5 



a 



+ 6) + (a - 6) sin 2 



2 iM 2 



(a + 6) — (a — 6) sin 2 0' 

and taking the square root of both sides gives us (13.1 5ft . □ 
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Now we can complete the proof of Gauss' theorem. We take the differential of the left 
hand side of (13.91) : we obtain 



, 2 cos 0'v/a 2 cos 2 0' + 6? sin 
cos (bad) = 



2 JJ 



a + b+ (a-b) sin 2 0' 

where we applied (|3.14p . 

Taking the differential of the right side of (13 .9p we get 



d 



2a sin 0' 1 2a cos 0'{(a + b) — (a — b) sin 0'} , 



a + 6 + (a - 6) sin 2 0' J {(a + 6) + (a - 6) sin 2 0'} 2 

Equating the right hand side of the previous two equations we and using (I3.15P we obtain 



2 cos 4>'^/a\ cos 2 0' + 5 2 sin 2 0' 2acos0'{(a + 6) - (a - b) sin 2 0'} 



a + 6+ (a -6) sin 2 0' {(a + 6) + (a - 6) sin 2 0'} 2 

JJ. a , IA / l\ ■ 2 -n {(a+fe)~(a-b)sin 2 ^} 2 

# a{(a + 6) - (a - 6) sin } , a {(a+b)-(a-b) sin 2 <A'} 



^a 2 cos 2 + 6 2 sin 2 {(a + 6) + (a - 6) sin 2 0'} 2 ^ a \ cos 2 0' + &f sin 



2 J./ 



\/ a 2 cos 2 0' + 6f sin 2 0' 

This completes the proof of ( 13. 12ft . therefore of ( 13.131) . and therefore, of Gauss' theorem. □ 
4 The Simple Pendulum 

First, we define the dynamical system. It is an idealization of a real pendulum. 

Definition 4.1. The simple pendulum consists of a particle which is constrained to move 
without friction on the circumference of a vertical circle and which is acted upon only by 
gravity. We describe it mechanically as follows: 

• a massless inextensible rigid rod has a point-mass attached to one end; 

• the rod is suspended from a frictionless pivot; 

• when the point-mass is given an initial push perpendicular to the rod, it will swing 
back and forth in one vertical plane and with a constant amplitude; 

• there is no air resistance. 

The following properties of the simple pendulum are readily available in numerous text- 
books. For example, the standard work of Pars 

The nonlinear differential equation which models the motion of the simple pendulum is 

H + ?*.» = », (4.D 



where g is the acceleration due to gravity, / is the length of the pendulum, and 6{t) is the 
angular displacement, at time t, of the pendulum measured positively (counter-clockwise) 
from the vertical equilibrium position. 

The Period of the pendulum, T, is the time taken by a double oscillation, to and fro, 
and is given by the following famous formula 

T = 4,/l/ f - d l (4.2) 

]j 9 JO Jl- S i n 2 I gin 2 

where a is the maximum angular displacement of the pendulum. 

The formula (14. 2 p shows that the period is proportional to a product of a function of the 
length, /, alone, and the maximum angular amplitude, a, alone. That is, there is already a 
"separation of variables" in the formula for the period. 

The integral in (" 14 . 2 [) is a complete elliptic integral of the first kind and, as we already 
noted, cannot be evaluated by any finite combination of elementary functions. So we must 
find suitable approximative formulas for K. 

It is customary to write k := sin | so that the integral in (14.21) is 

z sin o 



\[^k 

The quantity k is called the modulus of K and | is called the modular angle. In our case, 
the modular angle is one-half of the maximum angular displacement of the pendulum, and 
we write (with an abuse of notation) K(k) = K(a). The complimentary modulus, k' ^ is 
defined by k 2 + k' 2 = 1 and the corresponding complete elliptic integral of the first kind is 
K(k') = K'. 

Expanding H4.3[) by the binomial theorem and integrating term by term, we obtain 

(1.1) 

z . zj. . ,j . . . z/t I 

n=l L J y 

This gives us the fundamental theorem: 

Theorem 4.2. The period of a simple pendulum of length I, oscillating through an angle 2a, 
is equal to 

T 



K(k) = 1 { 1 + V P" 3 - 5 ---^-!) ] V" X 
y ' 2\ ^ 2-4-6---2n J 




□ 
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The last formula in Theorem 14.21 comes from substituting the M AcLaurin expansion of 
sin ^ into the previous series, and rearranging in increasing powers of a. 

Thus, if a pendulum swinging through an angle of 2a makes iV beats a day, and if a is 
increased by 5a, then the formula N • j — 86400, where 86400 is the number of seconds in 
a day, shows that the pendulum will lose 



43200 



1 



l\K(a) K(a + 5a) 



beats a day. 

For example, a pendulum, which beats seconds when swinging through an angle of 6° 
will lose about ll| seconds a day if made to swing through 8°, and about 26| seconds a day 
if made to swing through 10° [6]. 

If we truncate the previous series expansions for the period, we obtain the following 
approximative formulas: 




The first formula, 




(4.5) 



T (4.6) 



is the Huygens formula or the small angle approximation for the period. It does not contain 
a and gives an approximation which is independent of the period. Indeed, it is the formula for 
the period of simple harmonic motion, (SHM), realized by a particle travelling in a circular 

path of with constant angular velocity 

Just how accurate is the Huygens formula? It seems worthwhile to cite the lower bound 
found by Thurston [T3] and the upper bound found by Pars [□] • 

Corollary 4.3 (Pars-Thurston) . The Huygens small-angle approximation satisfies 



Mn f 2ttJ 1 - V sn. n 



/orO^a^f. □ 
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As Thurston points out, if he were to use the Huygens formula to adjust the length of 
his grandfather clock, which has an amplitude of 5°, to beat seconds, the error bounds show 
that the clock would lose between 4 and 8 minutes per week. 

These same bounds show that the Huygens formula is accurate to within 1% of the true 
period T for a smaller than about 14°. 

The second formula, 

7fl +(5)"K)j <« 

tells us that in the correction for the amplitude of a swing, the period must be increased by 
the fraction | sin 2 ^ of itself. Thus, if a pendulum swinging through an angle of 2a makes iV 
beats a day, and if a is increased by 5a, the pendulum will lose approximately • sin a • 5a) 
beats per day PQ. 

The last formula _ 

(4.9) 





is due to Daniel Bernoulli and is, historically, the first published correction term [12] to 
the Huygens formula (14.61) . 



5 The AGM approximations to the period 

The previous section shows that the problem of finding an approximative formula for the 
period T of the simple pendulum reduces to the problem of approximating the complete 
elliptic integral K{k) = K(a). 

If we take a := 1 and b := k' = cos ^ in the formula (I3.2p for Gauss' integral I(a, b), we 
obtain 

/( 1, cos ^W(*)= r , d ± / f . # (5.1) 

V 21 Jo ^1 - S in 2 f sin 2 Jo v/l-PsinV 

Applying the AGM to this choice of a and b we obtain the following sequences: 

a = 1 b = cos - 

1 / a\ 2 a i ( a \ 5 

ai = - 1 + cos — = cos — Oi = cos — 

2V 2/ 4 V 2/ 

If, / a\5l 2 , a ( a\\ 

a 2 = - < 1 + cos — > o 2 = cos — cos — 

4\ V 2/ J 4V 2 

if 2 a ( a\\ 
= -|cos 4 + (cos -J 

n 2 -.z' in i i 

a / a\2 , 1 / a\2 / «\2 / a\s 



a 3 = 71 cos — + ( cos — ) ) b 3 = - <( 1 + ( cos - J > ycos -j [cos - 
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We can use either — or — as an approximation to — . 

In order to discuss the accuracy of these approximations, we recall some definitions from 
numerical analysis. See Hildebrand |7J. Each digit of a number, except zero, which serves 
only to fix the position of the decimal point is called a significant digit or a significant figure 
of that number. 

Definition 5.1. (a) If any approximation iV to a number N has the property that both 
N and N round to the same set of significant figures, and if n is the LARGEST integer 
for which this statement is true, then N is said to approximate N to n significant 
digits. 



(b) 



, , . true value — approximate value E(N) , 

R(N) = relative error := ^— = —-^ , (5.2) 

v ' true value N K ' 



where E = E(N) is the absolute error. 

The importance of the relative error is shown in the following result. 
Proposition 5.2. N approximates N to n significant digits if and only if 



< {|- (5-3) 



Proposition 5.3. If R = R(N) and 



R = R(N) := ^1 , (5.4) 

then _ 

p p 

R = and R = = . (5.5) 

l-R l+R V ' 

Now we are ready to present Ingham's results. 
Theorem 5.4 (INGHAM). Let R n be the relative error in the approximation — m — and r n 
be the relative error in the approximation — f=s — taken positively. That is, let 



1 \ 1 1 / 1 \ 1 



1 + r n J b n jj, \ 1 - Rn ) a. 



(5.6) 



Then, 



0<r n <R n < ^ • (5.7) 

Za n +i 
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Proof. If we divide the numerator and denominator of (12.111) by fi, we see that the relative 
error, R n , in the approximation fa a n is greater than the relative error, r n , in the approxi- 
mation n?ab n . Substituting the formulas (15. 5p into the inequality < r n < R n , some simple 
algebra leads us to the inequality 

< r n < Rn. 



Moreover, 



Rn 



- ~ — - -h 



1 



< 



Qj n b n Q, n b n 
o-n + b n+ i 2a n+ i 



(5.8) 
(5.9) 



where the second equality follows from 



2 j.2 



a n (a 



l n u n+l 

and the second inequality follows from 

a-n + b n+ i > a n + b n = 2a n+ \. 



□ 



Theorem 5.5 (Ingham). If To denotes the Huygens small-angle approximation to the true 
period, T , then, for < a < n, 



H 


f 2 ] 






1 






1 + (C0S|)5_ 






cos -(cos 


(l + r 2 ) 



where 



1 f a a\ 4 

< r 2 < ito < — j - sin — tan — 

2 6 cos a V 4 4 / 



Proof. If we take n = 2 in the error estimate (15.7ft . and using 

8a n+ i(a n — &„) = (a„_i — fr n -i) 2 ) 

we obtain 

a 2 -6 2 _ 1 (a!-6!) 2 _ 1 1 (a ~b ) A _{l(a-b)Y 



(5.10) 



(5.11; 



(5.12) 



< r 2 < R 2 < 



2a, 



2a3 8a 3 



2a 3 8a 3 (8a 2 ) 2 



2 6 a^ 



i 3 u, 2 



(sin f ) 8 < (sin f ) 8 



(sin f ) 8 



1^a\al 2% 4 2 2 6 (cosf) 4 cos. 



. a a 
. sin — tan — 

2 6 cos s V 4 4 



□ 



Corollary 5.6 (Ingham). IfO < a ^ f ; i/ien 



< r 2 < i? 2 < 



70000 

and i/ras tae approximation is correct to 5 significant figures in the worst case, a = ~. 



(5.13) 
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Proof. If we take a = ^ in (15. lip and note that 



and therefore 



a 7rV2-\/2 a n ^2 + ^2 

sin — = sin — = , cos — = cos — = 

4 8 2 4 8 2 



cos 2 f sin 2 | 1 



V2 + 1 y/2-1 2^2' 



we conclude 



I f . a a\ 4 y/2 1 
sm — tan — 



V2 



2 6 cos|V 4 4 / 2 9 (v^+l)6 2 9 (99 + 70\/2) 

99 - 70^ 



2 8 f 



{since (99) 2 - (70 V2) 2 = 1} 



2 8 2 . 70 • 2.000050 . . . 
1 1 1 

< ^ — = — — < 



2 10 ■ 70 71680 70000 ' 



(5.14) 
(5.15) 



(I) (I) 

Finally, we have shown that i?2 < -^r < —^r and this means that the approximation is 



correct to 5 significant digits. 



10 5 10 5 



□ 



The elegant calculations in this proof are due to Ingham [TT]. If we apply the above 
computations to the case n = 3, we obtain the following results. 

Theorem 5.7 (Ingham). IfT denotes the Huygens small-angle approximation to the true 
period, T, then, for < a < n, 













f 2 1 


r 


(l-ife) 


i-\ 


_cos 2 f + (cosf)i - 







and 



where 




< r 3 < R 3 < 



1 / . a a 
—, r— sm — tan — 

2 14 cos 2 s V 4 4 



Corollary 5.8 (Ingham). J/0 < ck < | ; t/ien 

< r 3 < i? 3 < 



1 



(5.16) 



(5.17) 



(5.18) 



20000000000 

and thus the approximation is correct to 10 significant figures in the worst case, a = |. 



(5.19) 
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If we take a = 179° in f )5.18p . we obtain R 3 < 4.50 ... % which is in fair agreement with 
the machine calculation of [2] whose machine calculations showed R3 ~ 1%. In the next 
section we show how to bring it into much closer agreement. 

6 Renormalization 

Let k := sin 75° = and let the complimentary modulus k' := sin 15° = v/ ^~ v ^ . Let 

K y K' be the corresponding complete elliptic integrals of the first kind. In 1811, A.-M. 
LEGENDRE proved the following remarkable result [S] : 

K = \fl-K'. (6.1) 

This equation implies that a pendulum with an amplitude of 300° and length / has the same 
period as a pendulum with an amplitude of 60° and a length 3/. 

Towards the end of the nineteenth century Greenhill [6] proved: let the center of the 
circle of the pendulum's trajectory be O. Let B'B and b'b be two horizontal parallel chords 
of length where B'B is above the center and b'b is below the center. Let k := sin -.ZB OB 

Id 4 

and k' := sin ^Zb'Ob. Then 

K = V7-K'. (6.2) 

This equation implies that a pendulum of length / and amplitude ZB'OB has the same 
period as a pendulum of length 71 and amplitude Zb'Ob. 

That is, by suitably decreasing the amplitude and simultaneously increasing the length, 
one obtains a new pendulum with the same period. The smaller the amplitude, the more 
exact is each approximative formula we have developed. Thus, if we perform this process 
of replacing a given pendulum with pendulums of longer lengths and smaller amplitudes, our 
approximations become better and better. 

This is an example of what today is called renormalization and it is a fundamental 
technique in the study of dynamical systems. The two examples, above, are taken from 
the deep and beautiful theory of complex multiplication of elliptic functions (curves), which 
is one of the most active branches of research on the frontiers of modern mathematics. 
Unfortunately, it lies outside the scope of our presentation (see [I] and [15]). 

It turns out that the AGM furnishes us with another example of pendulum renormaliza- 
tion. Let us look at the equation 

I(a,b) = J(ai,&i) (6.3) 
for the case a := 1, b := cos | which gives us ai = cos 2 j , b\ — cos |)^. We note that 

1 1 1 _ 1 1 

Va 2 cosV + & 2 sinV ~ «i ^/l _ | s in 2 0' ~ cos2 f yjl - tan 4 f sin 2 0' ' 
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Returning to the formula for the period T of oscillation in an angle of 4a, we obtain the 
equation 



(cos 4 f ) 



9J ° ; /l-sin 2 f sin 2 V 9 Jo ; /l - (tan 2 f ) 2 sin 2 0' 



9 J o Jl - sin 2 ^sinV 



This equation can be stated as follows. 



Theorem 6.1. One iteration of the AGM transforms the pendulum of length I and maximum 
angular displacement a into another pendulum with the same period but with new length 



cos 4 f 



(6.4) 



which is longer than the original length I, and new maximum angular displacement 

(6.5) 



a% := 2 • arcsin(tan 2 — ) 



which is smaller than the original. □ 

This theorem allows us to "explain" the results of the numerical investigations of the 
accuracy of the AGM which Carvalhaes and Suppes presented in [2]. One uses the 
AGM to renormalize (reduce) the angular displacement a so that the Ingham estimates are 
applicable. 

1 T 

Example 6.2. Carvalhaes and Suppes state that — approximates — to within 1% for 
< a ^ 163.10° while the Ingham bound ( 15. lip gives ^ a ^ 162.5°. Again they report 
that the approximation has a relative error no bigger than for ^ a ^ 4.57° while the 
Ingham bound 05.1 lj) gives ^ a < 4.258°. 

Example 6.3. As another example, Carvalhaes and Suppes state that — approximates 

T 1 

— to within 1% for ^ a ^ 179.99°. The Ingham bound (I5.18P gives that — approximates 
T a 3 

T 

— to within 1% for ^ a ^ 177.98°. But, one application of the AGM reduces a = 179.99° 
To 

to a = 177.85° and now the Ingham bound (I5.18P shows that three more applications of 

1 T 

the AGM give an approximation to within 1%, that is, — approximates — to within 1% 

a 4 T 

for ^ a ^ 179.99°, in agreement with [2]. 
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The remaining results in [2] can be "explained" similarly. 

We cannot emphasize strongly enough the importance of rigorous upper and lower bounds 
for the absolute and relative errors in approximative formulas. Our analysis allows us to 
predict "a priori" the accuracy of a given approximative formula as well as to justify the 
resulting numerical studies. Moreover, Ingham's elegant and beautiful investigations give us 
practical tools to tailor our approximative formulas to the needs of the accuracy demanded. 
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